Autocatalytic plume pinch-off 
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A localized source of buoyancy flux in a non-reactive fluid medium creates a plume. The flux can 
be provided by either heat, a compositional difference between the fluid comprising the plume and its 
surroundings, or a combination of both. For autocatalytic plumes produced by the iodate-arsenous 
acid reaction, however, buoyancy is produced along the entire reacting interface between the plume 
and its surroundings. Buoyancy production at the moving interface drives fluid motion, which in 
turn generates flow that advects the reaction front. As a consequence of this interplay between 
fluid flow and chemical reaction, autocatalytic plumes exhibit a rich dynamics during their ascent 
through the reactant medium. One of the more interesting dynamical features is the production of 
an accelerating vortical plume head that in certain cases pinches-off and detaches from the upwelling 
conduit. After pinch-off, a new plume head forms in the conduit below, and this can lead to multiple 
generations of plume heads for a single plume initiation. We investigated the pinch-off process using 
both experimentation and simulation. Experiments were performed using various concentrations of 
glycerol, in which it was found that repeated pinch-off occurs exclusively in a specific concentration 
range. Autocatalytic plume simulations revealed that pinch-off is triggered by the appearance of 
accelerating flow in the plume conduit. 



PACS numbers: 47.20.Bp, 47.70.Fw, 47.15.-x 

Plumes are organized flow structures of considerable 
importance [TJ [5] due to their widespread natural occur- 
rence, particularly in geophysical flows [3H6]. A plume 
can be steady, in the case when the flow structure is 
well-developed and is no longer time-dependent, or it can 
be starting, in the transient case where buoyant fluid is 
penetrating the medium into which it is growing. The 
majority of work on plumes has focused on the steady 
variety, since the persistence of the columnar conduit fol- 
lowing behind the head plays a significant role in many 
geophysical processes. Long-lived plume conduits are be- 
lieved to underlie the formation of mantle hot spots, such 
as the one responsible for the Hawaiian Island - Emperor 
Seamount chain [3] . We will consider autocatalytic start- 
ing plumes in this paper. 

Before a laminar starting plume reaches a steady state, 
it develops its most prominent feature: an evolving head 
encompassing a vortex ring. Usually, the plume head re- 
mains attached to the conduit and the conduit remains 
linked to the buoyancy source [7] . Vortex rings are the 
subject of a vast literature and have a rich history [5]. 
They can be produced experimentally by using a piston 
to transiently inject a finite volume of neutrally buoyant 
fluid into quiescent surroundings [3] , a technique is simi- 
lar to the vortex ring generator pioneered by P. G. Tait in 
the nineteenth century |10j . Pinch- off refers to the pro- 
cess by which the vortex ring detaches from its source. 
In addition to neutrally buoyant scenarios, vortex ring 
formation and pinch-off has also been investigated in the 
context of buoyant starting plumes, where results suggest 
that dimensionless circulation of the ring is a universal 
quantity, regardless of whether it is produced by pinch-off 
from a plume head or formed by a piston pTJ [12] . 

Essentially free buoyant vortex rings can also be cre- 
ated by the pinch-off of plume heads driven by an au- 



tocatalytic chemical reaction [T3]. Nonlinear chemical 
kinetics can produce fronts of reaction which are ad- 
vected by the buoyancy-driven flow created by the re- 
action itself jUj. This "self-stirring" can, in turn, in- 
crease the production of the product species and change 
the morphology of the reaction front. In a container 
with a relatively unconfined geometry, the localized ini- 
tiation of a nonlinear autocatalytic reaction creates a 
special sort of reacting plume [121 [TS]. Such plumes 
can be viewed as weakly driven versions of expanding 
flame fronts, such as the flame "bubbles" that detonate 
Type IA supernovae [16j . Unlike supernovae, autocat- 
alytic plumes are readily produced and studied in the 
laboratory [HI Q22 H7] . 

This paper describes experimental and numerical in- 
vestigations on the behavior of autocatalytic plumes pro- 
duced by the propagating front of the iodate-arsenous 
acid (IAA) reaction [HI [19], focusing on the pinch-off 
process. The IAA reaction produces buoyancy-driven 
flow by its exothermicity and by the isothermal density 
difference between its reactant and product solutions [2U] • 
We explore the effect of adding glycerol to the IAA solu- 
tion on autocatalytic plume morphology and dynamics. 
Addition of glycerol increases the viscosity of the solu- 
tion, and it slows the rate of diffusion of the autocata- 
lyst. We found that pinch-off behavior is most prevalent 
in a specific window of glycerol concentration. Plumes 
produced with glycerol concentrations outside this win- 
dow did not produce multiple pinch-off events. To com- 
plement the experiments, we constructed an axisymmet- 
ric numerical simulation of autocatalytic plume dynamics 
and used it to examine the development of the reaction- 
driven flow that leads to pinch-off. We conclude by iden- 
tifying the subtle processes in the head and conduit that 
lead up to pinch-off events. 
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FIG. 1: A schematic of the apparatus. The reactant solution 
is contained in one of two different cylinders used for the ex- 
periments. These cylinders had lengths of 32 and 40.6 cm, 
and both had inner diameters of 8.9 cm. The capillary tube 
had an inner diameter of 2.7 mm. 



I. EXPERIMENT 



A. Apparatus 



The apparatus used for producing IAA starting plumes 
is shown in Fig. [I] The reaction tank was a large glass 
cylinder sealed at each end by large rubber bungs. A 
capillary tube entered through a hole in the lower bung. 
The outside end of the capillary tube was sealed with 
a short plastic tube clamped at one end which formed 
the initiation volume. The plastic tube was filled by a 
porous plug made of loosely packed cotton. Reactions 
were initiated by inserting a thin hypodermic needle into 
the plastic tube and then injecting a very small amount 
of product solution into the porous plug. The purpose of 
the plug was to to quench any potential hydrodynamic 
disturbance of the quiescent reactant solution during re- 
action initiation. As a reaction front ascends, it must 
first make its way through the plug, and then into the 
capillary tube. Upon reaching the end of the submerged 
portion of the capillary tube, a reaction front is then 
able to escape into the large volume of reactant solution, 
launching a plume. The apparatus was illuminated from 
behind and still images of the evolving plume were cap- 
tured using a digital camera. 
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TABLE I: Fluid properties at 24.0°C, and the Schmidt num- 
ber, Sc, of the reactant fluids for various glycerol-water con- 
centrations, given in both volume percent %v and mass per- 
cent %m. D is the diffusion constant for the iodide autocat- 
alyst. 



B. Reactant preparation and fluid properties 

Reactant solutions were prepared using distilled water 
and reagent grade chemicals. Iodate stock solutions were 
made by dissolving KIO3 powder in distilled water. Ar- 
senous acid stock solutions were prepared by dissolving 
AS2O3 powder in stirred water heated on a hot plate. 
These stock solutions were diluted so that the reactant 
solution contained [IO3HO.OO5M and [As(III)]=0.020M. 
Congo Red indicator was added to all reactant solutions 
in order to visualize the reaction fronts. It was present 
in solution at a concentration of 2 x 10 _5 M. Congo Red 
changes color over a pH range of 3.0 to 5.0, where the 
acid form is blue and the base form is red. The reac- 
tion front leaves in its wake a product solution with pH 
of r-j 2.7, so that the upwelling blue products are eas- 
ily visible within the red unreaeted solution. In addition 
to the reactants, most solutions also contained chemi- 
cally inert glycerol, the effect of which was to vary the 
fluid properties of the reactant mixtures, which are given 
in Table [I] Unless otherwise stated, glycerol concentra- 
tions throughout this paper are given in volume percent. 
Densities, p, were measured with an Anton-Paar densit- 
ometer, and kinematic viscosities, v, were determined by 
interpolating data from Ref. [21] ■ Values for the diffusion 
constant D of the iodide I - autocatalyst were calculated 
using the Stokes-Einstein relation [22], starting from the 
measured viscosity and the known diffusion constant of 
I - in water [23] • All experiments took place at a room 
temperature of (24.0 ± 1.0)° C. 



C. Starting plume morphology and pinch-off 

After escaping the capillary tube, the advancing auto- 
catalytic reaction front generates a starting plume that 
accelerates upward and may exhibit one or more pinch-off 
events, depending on the experimental conditions. The 
various stages in the evolution of a plume that pinches-off 
were first described in Ref. |13) . In summary, these stages 




FIG. 2: (Color online) The morphology of autocatalytic chemical starting plumes in a) water, and in glycerol concentrations 
of b) 20%, c) 30%, d) 40%, and e) 50%. All images are of plumes at a height of ~20 cm from the outlet. It took the plumes 
a) 255 s , b) 555 s, c) 425 s, d) 449 s, and e) 681 s to ascend to this height before these images were captured. 



are as follows. In the initial stage the plume creeps slowly 
out of the capillary tube and its head remains roughly 
spherical. Following this stage, the plume head changes 
shape and the vortical motion within the head becomes 
visible as the product solution wraps within an overturn- 
ing vortex ring. As the plume head continues to ascend 
and enlarge, it begins to pinch-off from the upwelling 
conduit. During this stage of growth, the head becomes 
almost completely detached from the conduit, and be- 
comes an essentially free vortex ring. Beneath this vortex 
ring, a rising bulge in the conduit forms which eventu- 
ally becomes a second generation plume head. Multiple 
generations of plume heads can be formed through this 
pinch-off process. 

Whether or not an autocatalytic plume head pinches- 
off, and also the number of pinch-offs, depends on the 
fluid properties of the reactant solution. Fig. [5] shows 
the morphologies of fully developed starting plumes pro- 
duced in solutions with various concentrations of glyc- 
erol. Plumes produced in the solutions without glycerol, 
such as the one shown in Fig. [2^, did not exhibit pinch- 
off. Similarly, in 20% glycerol solutions pinch-off was not 
observed. All plumes that did exhibit pinch-off, such as 
those shown in Fig.[2]>e, generally had similar morpholo- 
gies, regardless of the number of times that pinch-off oc- 
curred. For experiments in 30% and 40% glycerol, there 
were always two or more pinch-offs. At glycerol concen- 
trations of 50%, however, pinch-off was rarely observed. 
In cases where it did happen, such as the run shown in 
Fig. [2]3, only a single pinch-off event took place. 

In addition to fluid properties, the number of pinch- 
offs observed for a given reactant solution composition 



also depended on the length of the cylinder in which ex- 
periments were performed. In the 32 cm high cylinder, 
the maximum number of fully developed plume heads 
produced was four. This occurred for a run in a 30% 
glycerol solution. To test whether or not the number of 
pinch-offs could be increased by giving plumes more ver- 
tical space to develop, experimental runs were performed 
in an apparatus with the same diameter as the 32 cm 
long cylinder, but with a length of 40.6 cm. These runs 
were only performed for glycerol concentrations of 20%, 
40%, and 50%. For 50% glycerol plumes, no additional 
plume heads were generated. However, for 20% glycerol 
plumes, extending the length allowed for a single pinch- 
off to occur, when, in contrast, no pinch-offs were ob- 
served in the shorter cylinder for the same solution. The 
maximum number of fully developed heads for a 40% 
glycerol run in the longer cylinder was five, compared to 
just three heads in the shorter cylinder. The ascent data 
for five plume heads produced in a 40% glycerol solution 
is shown in Fig. [3] 



D. Plume head Reynolds number 

For each experimental run, the ascent velocity of each 
plume head was calculated using linear fits to the head 
height as a function of time. The maximum head velocity 
Vmax was determined from the fit to the five consecutive 
data points giving the largest slope. From v max , the head 
Reynolds number, 
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FIG. 3: (Color online) The height as a function of time for five 
plume heads produced by a 40% glycerol run in the 40.6 cm 
long cylinder. Lines indicate linear fits to the five consecutive 
data points with the maximum slope. Slopes of lines such as 
these were used to determine v max for each plume head. 



was calculated, where Wh is the width of the plume head 
at the average time over which v max was measured. Wh 
was measured from digital images, taking into account 
the optical magnification due to the cylindrical shape of 
the reaction tank. We observed that Ret had similar 
values for all of the plume heads produced in a given 
run. The exception to this generality was the value of 
Re?! for the final head produced in a run with multiple 
heads, which was always somewhat less than the others. 
This was mainly due to the smaller value of its v max . 
Typically, for a given experimental run, the final plume 
head would develop close to the top of the tank, giving it 
less time to accelerate before it hit the top surface of the 
fluid. This limitation is enhanced by the growing mass 
of product solution at the top of the tank which further 
decreases the vertical space into which the final head can 
ascend. For all of our experiments, Reh never exceeded 
~ 20. This value was reached in a 20% glycerol run, 
where Wh ~ 2.9 cm, v rnax ~ 0.12 cm/s, and the kinematic 
viscosity of the solution was v = 1.75 x 10 -2 cm 2 /s. The 
value of Re^ indicates that while the flow around the 
plume head is certainly laminar, Re^ is not so small that 
it can be considered Stokes or "creeping" flow, which 
would require Rey t -C 1. 

To understand the general tendency of the plumes to 
pinch-off, we examined the trend in the number of pinch- 
offs and head Reynolds numbers for varying glycerol con- 
centration. The glycerol concentration is conveniently 
expressed in terms of the dimensionless Schmidt num- 
ber, 

Sc = v/D. (2) 

Recall that the glycerol concentration effects both the 
kinematic viscosity v and the molecular diffusivity of the 
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FIG. 4: (Color online) The average plume head Reynolds 
number (Reh) for glycerol solutions characterized by vari- 
ous Schmidt numbers, Sc. Diamond symbols indicate runs 
that did not pinch-off in the short cylinder apparatus. Circu- 
lar points indicate that pinch-off occurred, with shaded (red) 
points indicating conditions under which two or more pinch- 
offs were observed. 



iodide autocatalyst D. The complete dynamical equa- 
tions and their nondimensionalization are discussed in 
Section III Al below. 

For all experimental runs, the mean of the plume head 
Reynolds number Re?! for each run was calculated by 
averaging the Re^ values for all the heads, omitting the 
final one in the case of multiple pinch-offs. As discussed 
above, the last head often did not have the time to fully 
develop before interacting with the top of the tank. Re^ 
values from different runs were then averaged to obtain a 
characteristic Reynolds number (Re^) for all experiments 
that took place in fluids with same glycerol concentration 
and hence the same value of Sc. These data are shown 
in Fig. [4] Examples of the plume morphologies for each 
of the Schmidt numbers displayed in Fig. [4] are shown in 
Fig.0 

For the two smallest Schmidt numbers, corresponding 
to runs in 0% and 20% glycerol, (Reh) is approximately 
independent of Sc and no pinch-offs were observed. This 
reflects the rather small and undeveloped head morphol- 
ogy under these conditions. As Fig. [2^i shows, in pure 
water, the head remains small and compact. In 20% glyc- 
erol, a vortex ring develops in the head only near the top 
of the tank. While these plumes rise faster than plumes 
in solutions containing more glycerol, (Re^) is limited 
by the small size of Wh- As glycerol concentration in- 
creases, Sc increases and (Re^) begins to drop. In this 
regime, multiple pinch-offs are observed. This behavior 
is accompanied by a decrease in (RejJ with viscosity with 
only small changes in Wh and v max . Finally, at the high- 
est value of Sc, only one pinch-off is observed while (Re^) 
continues its downward trend. Thus, we find that mul- 
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FIG. 5: (Color online) The aspect ratio w^/Ph for autocat- 
alytic plumes in a 30% (dark circles, blue online) and 40% 
(light circles, orange online) glycerol solution. The horizontal 
grey line indicates the constant value of Wh/£h = 1-24 ± 0.04 
for non-reacting forced plumes |7j. Here, t = corresponds 
approximately to when a vortex ring first forms in the head, 
not to the moment when the reaction front emerges from the 
capillary tube. 
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FIG. 6: (Color online) The radius of curvature R c of the crest 
of an autocatalytic plume head as a function of time. The two 
straight lines are fits to the data before and after pinch-off. 
The inset shows the crest of an autocatalytic plume head and 
the best fit circle that describes its shape. The head was 
imaged at a height of 16.7 cm above the outlet, 448 s after 
the front emerged from the capillary. 



tiple pinch-off behavior occurs in a special, intermediate 
range of Sc. 

Plume heads reach maximum velocity because as the 
head grows in width, the flow field begins to have stronger 
viscous interactions with the cylindrical wall of the tank. 
As the head reaches widths comparable to the diameter 
of the cylinder, vortical fluid motion in the head appears 
to quench because of wall effects. Eventually, the con- 
strained vortical motion is no longer able to effectively 
entrain fresh reactant solution into the head, slowing the 
growth of buoyancy generated by the reaction. In the 
intermediate range of Sc, the plume rises slowly enough 
that a complex head morphology and multiple pinch-off 
events have time to develop. For the highest value of Sc, 
however, the increased constraint from boundary effects 
is likely the reason why pinch-off is seldom observed. 



E. Evolution of head shape 

In addition to the head width Wh, the head length £h, 
was also measured from digital images. We defined ih 
as the distance from the top of the head to the base of 
the underturning lobe beneath it. The dimensionless as- 
pect ratio Wh/(h is a useful way to compare autocatalytic 
plume heads to buoyant plume heads produced without 
chemical reactions [7J. For laminar forced plumes, the 
growing plume heads are self-similar and this ratio is 
constant [7J with Wh/ih — 1.24 ±0.04. In contrast, auto- 
catalytic plume heads do not remain self-similar. Fig. [5] 
shows that the aspect ratio of the autocatalytic plume 



head increases with time. A steep increase in the rate of 
change of Wh / £h occurs at approximately the same time 
that a narrowing of the conduit beneath the head be- 
comes obvious, indicating that a pinch-off event is under 
way. 

Another interesting feature of autocatalytic plume 
head morphology is the spherical shape of the top of the 
head, which we refer to as the crest. The inset of Fig. [6] 
shows an image of a 40% glycerol plume head with a cir- 
cle fitted to the crest shape. Positions extracted from the 
image using edge detection were corrected for the opti- 
cal distortion of the cylindrical tank and then fit to a 
circle to extract the radius of curvature R c . This analy- 
sis is similar to the method employed in a classic study 
of the shape of bubbles rising in a cylinder by Davies 
and Taylor gj. The growth of R c with time for a 40% 
glycerol plume is shown in Fig. [6j We find that R c goes 
through two distinct periods of linear growth. This is 
in contrast to the growth of the head width Wh, which 
has been shown in previous work to have three distinct 
growth regimes |13j . While the transition from a small 
head without vortical motion to one with a vortex ring 
marks a change in growth rate of the width, it does not 
cause a change in the growth rate of the crest curvature 
R c . Instead, the only transition observed in the growth 
rate of R c occurs after the head pinches-off from the con- 
duit and becomes an almost-free vortex ring. A steep 
increase in growth rate of R c is a very clear and unam- 
biguous morphological indicator of pinch-off. 
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II. AUTOCATALYTIC PLUME SIMULATION 

In this section, we turn to a discussion of the axisym- 
metric numerical simulation which was used to diagnose 
the processes leading up to pinch-off in autocatalytic 
plumes. While it was computationally intractable to use 
our numerical model to simulate fully realistic plumes 
with multiple pinch-off events, it was possible to use our 
model to gain insight into the main features of the pro- 
cess. 



A. Dimensionless scaling considerations 

To describe our numerical model, it is helpful to first 
focus on the various dimensionless groups which control 
the fluid mechanical, thermal, and chemical aspects of 
plume dynamics. Plume motion depends on three differ- 
ent diffusivities. The most rapidly diffusing quantity is 
the momentum, which scales with the kinematic viscos- 
ity v. The reaction front is exothermic, and the diffusion 
of heat scales with the thermal diffusivity n. Finally, the 
concentration of the autocatalyst is the slowest diffusing 
quantity, scaling with the diffusivity D. 

In addition to the Schmidt number Sc, defined in 
Eqn.[2j a second independent dimensionless ratio is the 
Lewis number, 

Le = k/D. (3) 

The Prandtl number, given by 

Pr = Sc/Lc = u/k, (4) 

is not required in our formulation, but naturally appears 
in discussions of purely thermal plumes pQ. 

For a 40% glycerol solution, using k = 1.2 x 10 -3 cm 2 /s 
and values from Table [TJ we find 

Sc = 9000, Le = 280, and Pr = 32. (5) 

The large values of Sc and Le indicate that momentum 
and heat diffuse much faster than chemical concentration, 
and hence the concentration profile at the reaction front 
is extremely sharp compared to all the other fields that 
surround an autocatalytic plume. In general, the wide 
range of Sc, Le and Pr indicate strong scale separation 
between the various effects. 

The IAA reaction produces buoyancy by means of both 
temperature and concentration changes. These are both 
sufficiently small that we may employ the Boussinesq ap- 
proximation, with the density of solution 

p = po[l - a T {T\ - T ) - a c (ci - c )], (6) 

where po is the initial density at reference temperature 
To and product concentration c\ = Co, with Co being the 
initial concentration of reaction product in the solution. 
a.T and a c are the thermal and compositional expansion 



coefficients, respectively. In this approximation, Eqn. [6] 
implies that the density is a linear function of the tem- 
perature and product concentration of the solution. 

Across the reaction front, there are both tempera- 
ture and concentration changes, which we denote by 
AT = T — T) and Ac = c\ — Co, respectively. In order to 
quantify the thermal and compositional contributions to 
buoyancy separately, we define two Rayleigh-like num- 
bers [551 HI]- The thermal and concentration Rayleigh 
numbers are 

Ra T = ^T i ! AT and R&c = ^ AC; (7) 
vD vD 

respectively. Here, L is a length scale that we choose 
to be equal to the diffusive thickness £ of the reaction 
front, defined by £ — \Dt, where r is a reaction time 
scale. We calculated £ by measuring the front velocity, 
Vf, in a capillary tube where convective effects were ab- 
sent, and by defining vj — \JDJt and r = i/vj. Pa- 
rameters for the simulation were for a 40% glycerol so- 
lution, in which we measured Vf = 9.9 x 10~ 4 cm/s. 
Using TJ = 4.2 x 10~ 10 m 2 /s from Table [TJ we found 
r = 4.3 s and I — 4.2 x 10~ 3 cm. These values, along 
with a T = 4.3 x 10~ 4 K" 1 and AT = 0.5 K , estimated 
from values given in Refs. [211 [27], gi ye R- a T = 0.10. We 
measured a c Ac = 3.6 x 10~ 4 , which along with the other 
relevant parameters gives Ra c = 0.17. Given the sim- 
ilarity between the values of Ra-r and Ra c , we see that 
the buoyancy effects from temperature and concentration 
changes are comparable in our system. 



B. Numerical model 

Our simulation of autocatalytic plume dynamics was 
based on a well-established model describing flow driven 
by an autocatalytic reaction in a thin vertical slot [231126] . 
The autocatalytic plume simulation was axisymmetric 
about the direction of gravity and carried out in the 
cylindrical coordinates (r, z) shown in Fig. [7j which also 
shows an example of the concentration field of a numer- 
ical plume. The full dimensionless equations for our au- 
tocatalytic plume model are 

V • u = 0, (8) 

i-^? = -Vj? + V 2 u + (Ra T T + Ra c c)e z , (9) 



-^=V 2 c + T(c), (10) 

and 

DT 

— =LeV 2 T + T(c). (11) 
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FIG. 7: (Color online) The cylindrical coordinate system 
used, with the concentration field for a plume produced by 
the simulation. 



Here, e z is an upward pointing unit vector and D/Dt is 
the material derivative. Equations |10| and 1 1 1 1 incorporate 
a nonlinear reaction term given by 
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FIG. 8: (Color online) The a) concentration, b) temperature, 
and c) vorticity fields for an autocatalytic plume that has 
pinched-off. The concentration profile shows a plume with 
c = 1 separated from the surrounding fluid with c = by a 
thin reaction front. The spatial domain for the simulation was 
rt = 500^ and Zj = 2000£, only data in the region bounded 
by r — 3001 and z — 1500 is shown. The time after initiation 
is 42. 5r. 



F(c) =c 2 (l-c), 



(12) 



Simulation results 



which describes the cubic autocatalytic IAA reaction. In 
this equation, the concentration c is dimensionless and 
comes from scaling ci, the concentration of the auto- 
catalytic product of the reaction, by its maximum value 
Cmax- The fully reacted fluid therefore has a dimension- 
less product concentration of c = C\j c max = 1, and com- 
pletely unreacted fluid has c = 0. The dimensionless tem- 
perature from the reaction T in eqs.[9]and 11 is bounded 
such that < T < 1, and comes from scaling AT by 
AT H = -AHAC / p Q C p . Here, AH is the heat of reac- 
tion and is negative for the exothermic reactions, and C p 
is the specific heat at constant pressure. 

The cylindrical computational domain was bounded by 
r = r&, 2 = and z = Zb- Non-slip boundary conditions 
were imposed on the velocity field, and no flow condi- 
tions on the energy and chemical concentrations. These 
conditions are given by 



dT dc 

u = — — = — = at z — and z = z^. (14) 
oz az 



The numerical calculation used the primitive velocity 
variables u = (u, v) and the non-hydrostatic pressure p, 
and was carried out using the "semi-implicit method for 
pressure linked equations" (SIMPLE) algorithm [2"B"I l2"Tj] . 



In this section, we compare our simulation results to 
the phenomenology of 40% glycerol plumes discussed pre- 
viously. We observed the growth and pinch-off of com- 
putational plumes, and all aspects of their dynamics and 
morphology were in good, albeit qualitative, agreement 
with experimental plumes. The main limitation of the 
numerical computation is in the details of the initial con- 
ditions and in the size of the computational domain, 
which was necessarily smaller than that of the experi- 
ment. An example of the concentration, temperature, 
and vorticity fields from the simulation of an autocat- 
alytic starting plume are shown in Fig. [8] 

The exact details of the inlet pipe used to launch 
plumes experimentally were not simulated. Instead, the 
initial conditions used for the simulation were a spheri- 
cal region of product with initial radius r = 15£ within 
which c = T=latt = 0. The spherical region was 
placed just above the bottom boundary so that its cen- 
tre was initially located at z = 16£ and r = 0. The full 
spatial domain used for the simulation had r& = 500£ and 
Zb = 2000£, and grid points with spacings of £. The time 
step was At = 0.005t. This computational domain cor- 
responds to a physical cylinder with a diameter of 4.2 cm 
and a height of 8.4 cm. These dimensions correspond to 
approximately one half of the diameter and one quarter 
of the height of the experimental apparatus. The com- 
putational expense was too great to extend the spatial 
domain beyond these dimensions. Although the region 
of the reacted fluid remained far from the boundaries, as 
it was for the experimental plumes, the region of nonzero 
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fluid velocity naturally extended to the boundaries and 
thus we may expect that the computational plume flow 
was more confined than in the experiments. 

The concentration field in Fig.|8^ shows a plume with a 
similar morphology to some of the experimental plumes 
shown in Fig. [2] At this stage of the evolution of the 
plume, it has already undergone the pinch-off process and 
a new, second generation head has started to form. The 
temperature field around the plume head in Fig. [8]d is 
much more diffuse than the concentration field compris- 
ing the head. This reflects the large value of the Lewis 
number, which is the ratio of these two diffusivities. In- 
terestingly, the warmest region of fluid is only partially 
enclosed by the reaction front, while most of it resides in 
the center of the vortex ring. This effect is similar to the 
localization of hotspots in the cusps of two dimensional 
fingered fronts [3D]. This central location also happens 
to be where the only positive vorticity appears. The vor- 
ticity field for a plume that has pinched-off is shown in 
Fig. [SJ;. There are two regions with negative minima in 
the vorticity field, located in the cores of the vortex rings, 
where there is strong clockwise fluid motion. The largest 
negative vorticity is in the vortex core of the first gener- 
ation head and the second minimum is located where the 
vortex ring in the second generation head is forming. 

Figure[9]shows that there is good qualitative agreement 
in the general behavior of experimental and simulated 
autocatalytic plume heads. Fig. shows that the simu- 
lated plume head accelerates as it ascends, even after it 
has pinched-off from the conduit, which was also observed 
experimentally. The pinched-off head continues to accel- 
erate until it reaches a maximum velocity, after which 
time it decelerates because of its interaction with the 
boundary at the top of the computational domain. The 
maximum velocity that the first generation head achieved 
in this simulation was v max = 67£/t = 0.065 cm/s. This 
compares reasonably well with v max — 0.090 cm/s ex- 
tracted from the experimental data for the first gener- 
ation plume head in Fig. [3]i. It is not surprising that 
the value of v max achieved by the simulated plume is less 
than that of the experimental plume because of the dif- 
ference in boundary constraints between the two. The 
simulated plume has greater viscous interaction with the 
boundaries, and it also has less vertical space to grow 
and accelerate. 

The simulation also shows good qualitative agreement 
between the morphological characteristics of simulated 
and experimental autocatalytic plume heads. Simulation 
results for Whjlh and R c are shown in Fig. [9]>c, which 
should be compared to Figs. [5] and [6] respectively. 
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time (r) 

FIG. 9: (Color online) a) The height of first and second gen- 
eration simulated plume heads as a function of time, b) The 
width to length Wh/ih aspect ratio and c) the radius of curva- 
ture R c of the crest of the first generation head. The shaded 
region in (c) shows the constant aspect ratio of non-reacting 
forced plumes [7J. Parts a), b), and c) show simulation data 
that correspond to the experimental data presented in Figs. 
[3j[5j and[6j respectively. 



D. Simulated pinch-off 

In their early stages, the external appearances of au- 
tocatalytic plumes are similar to non-reacting, laminar 
plumes [2 O [7] ■ This suggests that the interior flow pro- 
files of autocatalytic and non-reacting laminar plumes 



might be similar, at least initially. The simulation re- 
sults, however, show that this is not the case. A non- 
reacting laminar plume head is formed by buoyancy flux 
from its upwelling plume conduit. In order to feed the 
head, the conduit flow must rise at a greater velocity 
than the head [7J. In this scenario, the plume head can- 
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time (r) 

FIG. 10: (Color online) A space-time diagram showing the 
development of the center line velocity v c of a simulated plume. 
The units of v c are Ljr. 



not detach from the conduit. As the head rises, it must 
force its way up through quiescent fluid, and this induces 
large scale return flow in the tank. Once this flow has 
been established below the plume head, fluid rising in 
the conduit experiences less resistance to motion than 
was initially experienced by the head — the conduit is 
moving upwards into fluid that already has momentum in 
the direction of its motion. This leads to the formation of 
a steady-state conduit flow below the plume head, where 
the velocity profile across a horizontal cross section of the 
conduit is constant, independent of height |15j . 

The simulation results show that the transient phase 
during the development of the autocatalytic steady-state 
conduit is much more complex than in non-reacting lam- 
inar plumes. Previous experimental work [15 has shown 
that the conduit is conical, rather than cylindrical for the 
case of non-reacting laminar plumes. This shape is also 
clear in the simulated autocatalytic plume conduits. The 
most dramatic feature is of course the pinch-off process. 
In previous work, it was presumed that pinch-off occurs 
because of the acceleration of the autocatalytic plume 
head, which then leads to a "bottleneck" in the conduit 
at the site of pinch-off, resulting in the formation of a sec- 
ond generation head [13] . Our simulation results reveal 
a much more complete picture of the flow development 
involved in this process. 

As an autocatalytic plume evolves and the conduit 
elongates, reaction at the interface between reacted and 
unreacted fluids causes the buoyancy flux to be dis- 
tributed along the entire length of the conduit. The 
time-dependent flow in the conduit is shown in Fig. |10| 
which shows a space-time diagram of the vertical veloc- 



ity profile along the axis of symmetry of the plume. In 
the early stages of plume growth, before about 23r, the 
vertical velocity in the conduit increases with distance 
from the base of the plume. This shows that, even in 
the early stages, an autocatalytic plume conduit docs 
not play the same role as a non-reacting plume conduit, 
which delivers buoyancy flux to the more slowly ascend- 
ing non-reacting plume head. The ascent velocity of the 
autocatalytic plume head in the early regime exceeds the 
vertical velocity in the conduit, so that upwelling flow in 
the conduit does not reach the head. Instead of being 
fed buoyant fluid by the conduit, an autocatalytic plume 
head is able to generate its own buoyancy flux by means 
of reaction with the surrounding fluid. 

The later stages of autocatalytic plume development, 
leading to pinch-off, are more complex. Beyond approx- 
imately 23t, a velocity minimum becomes apparent be- 
neath the head. In Fig. |10[ the minimum is seen as a 
streak that trails below the maximum near the crest of 
the plume. The development of the velocity minimum 
occurs long before morphological indicators of pinch-off 
are observed, namely, the changes in the growth rates of 
either Wh/th or R c , or the visible narrowing of the con- 
duit beneath the head. All of these indicators suggest 
that pinch-off occurs from anywhere between approxi- 
mately 37r to 41t. The first generation head continues 
to accelerate upward, despite the development of the ve- 
locity minimum, highlighting the lack of influence that 
the conduit has on the subsequent evolution of the first 
generation autocatalytic plume head. The second gen- 
eration plume head is created by the accumulation of 
buoyant fluid rising up into the area of the velocity min- 
imum. After it is formed, the second head also acceler- 
ates away from the conduit, and a new velocity minimum 
develops in its wake, bringing the process of head gener- 
ation back to the beginning of the cycle. The simulation 
clearly shows the precursor structure in the conduit ve- 
locity that long precedes the morphological indicators of 
pinch-off. Although it superficially resembles the con- 
duit of a non-reacting plume, the structure below the 
head of an autocatalytic plume, which might better be 
called the "tail" , does not contribute in the same way to 
plume head motion. This is already the case long before 
pinch-off becomes apparent. 



III. CONCLUSIONS 

We have described an experimental and numerical 
study of buoyant laminar plumes driven by the auto- 
catalytic iodate-arsenous acid reaction. Experimentally, 
autocatalytic plumes were created in solutions where the 
viscosity was varied by the addition of glycerol. In the 
most viscous solutions, pinch-off of the plume head from 
the conduit was rarely observed, while in the lowest vis- 
cosity solutions, pinch-off was not observed. In the mid- 
range of glycerol concentrations used, however, pinch-off 
occurred multiple times from a single initiation event. 
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The successive heads all tended to achieve the same up- 
ward velocity. The morphology of autocatalytic plume 
heads was quantified in terms of their radius of curva- 
ture, R c and their width to length aspect ratio, Wh/£h- 
We found that this ratio evolved with time quite differ- 
ently from that of non-reacting forced plumes. A steep 
increase in the rate of change of R c was found to occur 
at the approximate time when the plume head visibly 
separates from the conduit. 

To further examine the pinch-off process, axisymmet- 
ric simulations of autocatalytic plumes were implemented 
for the fluid parameters corresponding to a reactant so- 
lution with 40% glycerol. Similar results were shown for 
both simulated and experimental plumes - first genera- 
tion plume heads were observed to accelerate and pinch- 
off from the conduit, and the metrics used to quantify 
the morphology of the first generation plume head were 
in good qualitative agreement with the ones measured 
for the experimental system. More importantly, the sim- 
ulation results provided new insight into the pinch-off 
process. It was found that subtle changes in the velocity 
structure in the conduit occur long before the appear- 
ance of external morphological changes associated with 



pinch-off. Furthermore, it was determined that the verti- 
cal component of the fluid velocity in the conduit is less 
than the ascent velocity of the head at all times, so that 
the conduit does not contribute buoyancy to the head but 
rather forms a tail in its wake. The appearance of veloc- 
ity minima in this tail are early precursors to the eventual 
pinch-off of the head. This velocity profile reveals that 
autocatalytic plume heads are nearly independent flow 
structures that are driven mainly by the chemical reac- 
tion between fluid in the heads and their surroundings. 
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